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Abstract 

Using the approach of low-energy effective field theory, the phase diagram is studied for a mixture 
of two species of pseudospin-^ Bose atoms with interspecies spin-exchange. There are four mean- 
field regimes on the parameter plane of g e and g z , where g e is the interspecies spin-exchange 
interaction strength, while g z is the difference between the interaction strength of interspecies 
scattering without spin-exchange of equal spins and that of unequal spins. Two regimes, with 
\dz\ > |<7e|) correspond to ground states with the total spins of the two species parallel or antiparallel 
along z direction, and the low energy excitations are equivalent to those of two-component spinless 
Bosons. The other two regimes, with \g e \ > \g z \, correspond to ground states with the total 
spins of the two species parallel or antiparallel on xy plane, and the low energy excitations are 
described by a sine-Gordon model coupled with a free scalar field, where the effective fields are 
combinations of the phases of the original four Boson fields. In (l+l)-dimension, they are described 
by Kosterlitz-Thouless renormalization group (RG) equations, and there are three sectors in the 
phase plane of a scaling dimension and a dimensionless parameter proportional to the strength of 
the cosine interaction, both depending on the densities. The gaps of these elementary excitations 
are experimental probes of the underlying many-body ground states. 
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I. INTRODUCTION 



Sine-Gordon model is important in field theory and statistical mechanics, and its renor- 
malization group equations define the Kosterlitz-Thouless universality class for a class of 
(l+l)-dimensional quantum systems and two-dimensional classical systems In this 

paper, we show that a sine-Gordon model coupled with a free scalar field emerges in the 
phase dynamics of a mixture of two different species of pseudospin-| Bose gases with in- 
terspecies spin-exchange interaction. Interestingly, here the scalar field described by the 
sine-Gordon model and the free scalar field are two different combinations of the phases of 
the four bosonic fields in this system. 

A mixture of two different species of pseudospin-| Bose gases with interspecies spin- 

spinor Bose gas as 



10] . In this model, 



exchange interaction exhibits novel features beyond a single species o 
well as a mixture of two species without interspecies spin exchange 
each atom has an internal degree of freedom represented as a pseudospin with basis states 
| t) an d | i), while there are two species of atoms, with the atom number of each species 
conserved. It can be described by the following Hamiltonian density, 

t 

era 

where a = a,b represents the two species and o =t, 4- V = V(x) is the external poten- 
tial, g , \ g } and g e are the interaction strengths for intraspecies scattering, interspecies 
scattering without spin exchange, and interspecies spin-exchange scattering respectively, 
proportional to the corresponding scattering lengths. For pseudospin- 1 atoms, intraspecies 



scattering strengths with and without spin-exchange are the same 11] . For simplicity, we 
assume g^? = g a for any a and a', gffl = g^f* = g s and g^ = g^ = g&. We define 
S a %{x) = ^aS^a, where ty a (%) = (jP^( x )i '4 ) i( x )) T , s * — r i/2, r* being the Pauli matrix, 
(z = x, y, z). Then TL can be rewritten as 



(2) 



~\~29e{S ax Sbx + SaySby) + 2>9z<Saz<Sbz, 

where g a b = 9s + gd, 9% = 9s — 9d- It can be seen that g a , and g a b characterize the usual 
density- density interactions, while g z and g e characterize the spin coupling between the two 
species. 



We make the presumption that g a > 0, g>& > and ^g a 9b > glpi which is needed for the 
stability of the system and can be naturally satisfied in reality |l0|. We study the phase 
diagram in the space of the parameters g e and g z , by using the approach of low energy 
effective field theory. The regime of g e > g z > has been discussed previously for higher 
dimensions, i.e. when the phase fluctuation is sup pre ssed such that the cosine of a phase 
variable can be approximated up to second order lOj. In this paper, we first extend the 
discussion to other parameter regimes. In the regime of g z > \g e \, the ground state is with 
the total spins of the two species antiparallel along z direction. In the regime of g z < —\g e \, 
the ground state is with the total spins of the two species parallel along z direction. In the 
regime of g e > \g z \, the ground state is with the total spins of the two species antiparallel 
on xy direction. In the regime of g e < — \g z \, the ground state is with the total spins of 
the two species parallel on xy direction. Then we focus on the case of \g e \ > \g z \ in (1+1)- 
dimension. Without approximating the cosine interaction term, the low energy excitations 
can be described by a sine-Gordon model coupled with a free scalar field, both fields being 
combinations of the phases of the original four boson field. It turns out that for given 
g e and g z with \g e \ > \g z \, there are three phases according to a scaling dimension and a 
dimensionless parameter proportional to \g e \. Both these two parameters depend on the 
densities of the two species. 



II. PHASE DIAGRAM ON g e - g z PARAMETER PLANE 

There is a symmetry between parameter points (g e , g z ) and (— g e , g z ). Consider the trans- 
formation i/;'^ = —i/Jat, ip'a\. = Vvu V'bt = = ^K- The Hamiltonian density in terms of 
the primmed operators in the parameter point (g e , g z ) has the same form as the Hamiltonian 
density in terms of the unprimed ones in the parameter point (—g e ,g z ). 

Now consider the case of g z > \g e \. Then from the Hamiltonian density (j2J), it is easy to 
see that in the ground state, S a and Sb must align oppositely in the z direction. We choose 
the mean field values in the ground state to be with ip° aJ[ = y/n^, ipai = 0> ^bt = 0> V'm, = \/^& 
so that S a = ^z and Sf, = —^z, where z is the unit vector in the z direction, n a is the total 
density of species a, (a = a,b). The low energy dynamics is dominated by the fluctuations 
of V'at an d il>bi, as the fluctuations of V><4 and ip^, whose mean field values are zero, must be 
of the amplitudes rather than the phases, and thus increase the energy. 
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Similarly, in the case of g z < —\g e \, the ground state is that with S a and <S& parallel in the 



z direction. We can choose the mean field values in the ground state to be with i\P a ^ 



ipai — 0, ^m- — s/nbi ^bi = 0j so ^at S a = ?j-z and S b = ^z. The low energy dynamics 
is dominated by the fluctuations of i/j^ and ipbt, as the fluctuations of ip a ^ and ipbi, whose 
mean field values are zero, must be of the amplitudes and thus increase the energy. 

In these two cases, which can be represented in a unified form as \g z \ > \g e \, the system 
behaves like a mixture of two species of spinless Boson gases, with the effective Hamiltonian 
density 



2m n 2 2 2 



whose excitation spectra are 



12| 



^ = 2 ( £ l + e D ± 2 V ~ £ D 2 + 4E a E b n a n b (g ab - \g z \Y, 



where we have introduced 



e 2 a = E a (2g a n a + E a 



(3) 



(4) 



(5) 



with a — a, b, and E a = being the free particle energy of species a. All spectra are 



2m a 



hat is, to — > when k — > 0. 



gapless, as in the usual case of phonon-like Goldstone modes. 

In the case of g e > \g z \, which has been discussed previously [l0|, the ground state is that 
with S a and S b antiparallel on the xy plane. One can choose the ground state to be with 

^ = ^ = ^j2, ^ = -vft = vW72- 

Similarly, in the case of g e < — \g z \, the ground state is that with S a and S b parallel on 
the xy plane. One can choose the ground state to be with 1/;^. = ip^ = ^Jn a /2, tp b ^ = tp b ^ = 

In the latter two cases, which can be represented in a unified form as \g e \ > \g z \, the 
effective Lagrangian describing the phase fluctuations is, with \g e \ replacing g e in the result 
for g e > g z > [lOj, that is, 



C 



eff = \(dtr T )A-\d t r) - ^(Vr^M-^Vr) + ^n a n b cos (2 74 ) 



(6) 



where 



72 

73 

V7 4 y 



4= -i- 
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(7) 
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with <3> a(T being the phase of ip acr , 

( 



A 



2 9a 9ab - \9e 



9ab ~ \9e\ 





2<?b 









\g e \V+ + 9z 

\9e\V- 






\9e\V- 
\9e\V+ ~ 9z J 



(8) 



with 



with 



2 n a n b 
( Jk o 
0^00 

o o e+ e- 
V o o £_ e+ y 



2 m a m 6 



(9) 



In (3+l)-dimension or (2+l)-dimension, the fluctuation of 74 is largely suppressed and 



we can make the approximation that 005(274) 



spectra can be obtained as 
k 2 



2jI, subsequently, the four excitation 



"hi 



9an a + g b n b 



m b 
1 



gmb_, 2 , (^6 - |#e|) 2 n a n 6 



+ 



m a m b 



IILIV 



(10) 
(11) 



+ A 2 =f VCA; 4 + DA; 2 + A 4 , 

where A 2 = \g 2 (^ + ^) - 2\g e \g z \n a n b , B = ^A(^- + *»-), C = 4(— ~ —? + g 2 z Jh ^, 
D = |c/ e |n a nJ^ 2 (^ - - 2*) - 2|^ e |^(^ + ^) + 2# 2 (^- + Under the conditions 

g a > 0, 4:g a g b > g 2 ah and |g e | > \g z \, all these excitations have real energies for any k, 
guaranteeing the stability of the ground state. 

It can be seen that uiy has a gap A while the other three excitations are gapless. That 
is, as k — » 0, Uijijn — > 0, but uiv — > A. 

Therefore in (3+l)-dimension, we obtain the mean-field phase diagram as shown in Fig. [IJ 
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FIG. 1: Mean-field phase diagram on g e — g z phase plane in (3+l)-dimension. In regime A, 
g z > \g e \, the ground state is V>° t = \A^> = 0, ip^ = 0, ip%± = y^. In regime B, g z < -\g e \, 
the ground state is tp®^ = y 7 ^, ip^i = 0, tjJL = y'rtb, ^{Ji = 0- I* 1 regime C, g e > \g z \, the ground 
state is ifi®*. = ip^> = yJn a /2, ip^L = —ip^ = \fn\j2. In regime D, g e < —\g z \, the ground state is 
€t = V'a ; = yfafi, < = < = V^2. 

III. RENORMALIZATION GROUP ANALYSIS IN THE CASE OF \g e \ > \g z \ IN 
(l+l)-DIMENSION 

Reconsider the case of \g e \ > \g z \. In (l+l)-dimension, the fluctuation is important and 
we must take into account the whole effect of cos(27 4 ) term, which is the only interaction 
term in £ e //> where 71 and 7 2 are both free fields and not coupled to 74. Hence we can focus 
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on the sine- Gordon field 74 coupled to a free scalar field 73 

<9 t 73 <%74 



C 



1 



1 



-4 ( ^73 ^74 





9e\V+ +9z 

+ -^-n a n b cos (274). 



(12) 



We now make a renormalization group (RG) analysis of the cosine interaction term due 
to spin exchange in (l+l)-dimension, by following the approach in jl|. It turns out that it 
still belongs to the Kosterlitz-Thouless universality class. But the novelty is that the scalar 
fields are now combinations of the phases of the original bosonic fields. 

Let us define x° = vt, with v = \ hnr^T^r^-, and two dimensionless field variables 
corresponding to 73 and 74 



X 



r £+(\9e\V+ + 9z) ,1 

[ %9l - 9l) 1 ' 
^+{\9e\v+ + 9 z),\ 

[ m-9i) 



l 74- 



Then the action can be written as 



S — Sq + Sj, 



(13) 



where So = J CocPx, Si = j Cid 2 x, with 



1 



d <p d x 




dup dxX 



1 p d\ip 



p 1 



dix 



Hi = ^cos {fix). 



(14) 
(15) 



where 



h = 

f2 = 



\9e\V+-9z 

\g e \v+ + 9z'' 





9e 


V- 




9c 


V+ + 9z 



V 



P = 2[ 



£+(\9e\V++9z) -,-i 



A = 



\g e \n a n b a z 
2v 
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a is the short range cut-off, which is the coherence or healing length, and can be estimated 
to be h[mb{g a n 2 a + g^nf + (g s + gd — ge)n a nb)}~ 1 ^ 2 , assuming m a > rrib- A is dimensionless. 
From £ the free propagator of x is obtained as 

Q (k] = fiK + fc| = 1_ Acos^ + sin^ 

x[ ' k 2 {hkl + k 2 )-{f 2 kl+pkl) 2 k 2 f 1 cos 2 9 + sm 2 9-(f 2 cos 2 9 + psm 2 9) 2, 

(16) 

where k 2 = k" 2 + kf, tan 6 = |jk 

Since the interaction term does not involve <p field, we only need to split \ into the fast 
and slow components, 

Xa(x) = Xa'(x) + h(x), (17) 



where 



Xa>(x) = £ e lkX Xk, (18) 



fc<A' 



Kx)= J2 eikX Xk, (19) 

A'<fc<A 

A = i is the momentum cut-off, A' = A — dA. The partition function can be written as |l| 

Z A = J VipVxA>Vhe- So[ *> x * ] - So[h] - Slbc * +h] (20) 
= Z h J V V V X A'e- So[m] (e- Sl[xA ' +h] )h, (21) 

where Z h = JVhe~ s °^, (...) h means taking average over the fast components of x- 
The effective action is thus 

S[<p, X A>] = S [^XA']-He- Sl[xA ' +h] )h 

~ S [<p,XA>) + (SiiXA> + h)) h - l((S][ X A' + h]) h - (S z \xa> + h]) 2 h ), (22) 

which allows us to calculate the RG flows of A and (3. We have 

(h(x)h(0)) h = f ^rM^ = ^T dl ' ( 23 ) 



where 



dl = T , 



K ( Ar ) = f *fL /iCos 2 g + sin 2 fl iArcos , 

1 ; "J 2n f lC os 2 9 + sin 2 9- {f 2 cos 2 9 + psin 2 9f ' 1 } 



d9_ f\ cos 2 9 + sin 2 9 

2^ f\cos 2 9 + sm 2 9 - ( 

which is a measure of correlation of fluctuations. 

8 



Then 

{eW'hn = e -^< fca W>* = 1 - £^dl = 1 - Ddl. (25) 
It is seen that the coupling between ip and x modifies the scaling dimension of e ll3x from 

which means the interaction term cos(/3x) is less relevant than that in the pure SG model, 
since k(0) > 1. 

Now we can obtain the renormalized action. Following [1], we have 

72 

(Si\xa> +h]) h = A(l - DdZ) | cos(/3 X aO, (27) 
(S?>* - (S,)* = -aA a £>gdi / d 2 x(V X A>) 2 , (28) 

where 

a = f^z 3 K(z). (29) 
Therefore, by rescaling A' — > A, the purely x-dependent part of the action is renormalized 

to 

S[ X ) = 1(1 + a\ 2 D 2 dl) J d 2 x{V X f + A[l + (2 - Ddl)) J ^ cos(/3 X ). (30) 

The overall factor in front of the Gaussian part of the action requires a renormalization of 
the field Xi as we U as the parameter /3 such that fix is invariant. That is, 

x{x) ^ (l + a\ 2 D 2 diy/ 2 X (x), 
-> (1 + a\ 2 D 2 dl)- l ' 2 p, 

which gives RG flows of A and j3. 

There is another field ip coupled to x- The coupling terms are those proportional to f2 
and p, respectively, in (114)) . The above renormalization of x an d lead to the renormaliztion 
of fi and p. As in the pure SG model, we introduce 



t = D - 2 



and assume t and A are small, as the parameter point of t = and A = is a fixed point of 
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the RG flows. Then to the order of A 2 , the RG equations read 



dt _ SaX 



dl k 2 (0) ' 

— = — \f 

dl ~ ' a (32) 

dl /?(0)^ 2 ' 
dp _ 2aX 2 

dl K 2 (0)^- 

IV. PHASE DIAGRAM ON THE t - y PARAMETER PLANE AND THE MASS 
GAP IN (l+l)-DIMENSION 

From the RG equations, one obtains 

d(t 2 ) - -^-A* 2 ) = 0, (33) 

which is similar to the equation in the pure SG model 1|, except that a and k are not 
constant here. Moreover, and % are both proportional to A 2 . Hence to the order of 
A 3 , one can replace ^f^y<i(A 2 ) as rf(^|p-A 2 ). Thus we arrive at the following equation, 

t 2 - y 2 = /x 2 , (34) 

where 



V = ~k{0) (35) 
/i represents a constant. For given g e and g z with \g e \ > \g z \, Equation ( 1341 determines the 

phase diagram of the model on the plane (t, y) in the regime where t and y are small, as 

schematically shown in Fig. [2j 

It can be seen that the (t, y) phase space is divided to three sectors, namely, weak coupling, 
strong coupling and crossover sectors. In the weak coupling sector, the effective theory scales 
to a Gaussian model, y(l) — > as I — > oo, and the spectrum is massless, while in the crossover 
and strong coupling sectors, the coupling constants flow away from the Gaussian fixed line 
and the spectrum has a mass gap. 

Note that both t and y depend not only on the interspecies spin-exchange coupling, but 
also on the densities. Consequently the phase is dependent not only on the interaction 
strengths, but also on the densities of the two species, which can be easily adjusted in 
experiments. To illustrate this explicitly, let us set n a = = n so that U = jjtg, f 2 = 0, 

1 3 

p = ^~^ a and thus k and a as well, are all independent of n, while /3 oc n"« and A oc n?. 
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o t 



FIG. 2: Phase diagram of our emergent sine-Gordon model coupled with a free field in (1+1)- 
dimension, with y = ^^A, t = D — 2, D is a scaling dimension. There are two separatrices t = ±y 
that divide the phase plane into three sectors: (1) t > y, the weak coupling (WC) sector; (2) 
\t\ < y, the crossover (C) sector; (3) t < —y, the strong coupling (SC) sector. The RG flows are 
similar to the pure sine-Gordon model. 

Therefore D oc and y oc ni. Then according to Fig. (2], for small enough n, the system is 
in the weak coupling phase. For large enough n, the system is in the strong coupling phase. 
Therefore, following the change of density, the system goes through phase transitions. 

The mass gap in the crossover and strong coupling sectors can be qualitatively obtained. 
In the strong coupling sector, fi is real, while in the crossover sector, [i is purely imaginary. 
Following 1], it can be found that the mass gap in the strong coupling sector is 



M 



Uo) (36) 



Aexp(-l/y ), V < \yo\, 

where the subscript "0" means the bare values, that is, the values measured in experiments, 
while in the crossover sector, the mass gap is 

M = Aexp(-7r/2|/i|), (37) 

with | /x | ^> \t \. 

The scaling behavior of the mass gap may be observed in experiments. Taking M = 



Aexp(— lA/o) as an example. If n a = n& = n, we have = ~ nz, then the relation 

between M and n may be investigated. 

The last two of the RG equations (l32il determine the RG flows of the couplings between 
the fields if and x- Since a > 0, it is easy to find that f 2 = 0, p = is the only stable 
fixed point of the two equations, namely, whatever the initial values of f 2 and p are, they 
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inevitably flow to 0. Moreover, the larger A is, the more rapidly f 2 and p flow to 0. It is like 
that its strong self interaction "traps" the field \ and separate it from (p. If the bare value 
of A were 0, there would be no RG flows of f 2 and p. 

Also note that if we diagonalize Cq in (|14p . then Cj in (|T5|) becomes a cosine term of 
cosine of a linear combination of the two fields, of which the RG analysis is quite difficult. 
Hence we use the above approach instead. 

The elementary excitations studied here can be experimentally measured by using the 
Bragg spectroscopy. The gap in a collective mode is a novel feature absent in the BEC 
mixtures previously studied. The two key parameters g e and g z both originate from the 
interspecies spin- dependent scattering, thus they are roughly of the same order of magnitude. 
We expect that the excitation gaps can be detected in experiments and are indications of 
the underlying many-body ground states. 

V. SUMMARY 

We have developed a low energy effective theory for a mixture of two species of pseudospin- 
| Bose gases and explore the phase transitions in the space of the parameters g e and g z , 
where g e is the interspecies spin-exchange interaction strength, while g z is the difference 
between the strengths of equal-spin and unequal-spin interspecies interaction without spin 
exchange. The phase diagram on the plane of parameters g e and g z is shown in Fig [TJ In 
the regime of \g z \ > \g e \, the system is effectively described by a two component model, and 
the excitation spectra are gapless. In the regime of \g e \ > \g z \, the system is described by a 
four effective fields, which are combinations of the phases of the four original boson fields. 
There is a cosine interaction term of one of the effective field, which can be approximated 
as a square in (3+l)-dimension. There are three gapless modes and one gapped mode. 

In (l+l)-dimension, the effective theory in the regime of \g e \ > \g z \ is a novel realiza- 
tion of a sine-Gordon model coupled with a free scalar field, on which we have made a 
renormalization analysis. Described by Kosterlitz-Thouless equations, the phase space is 
further divided into three sectors, as shown in Fig. [2J according to a scaling dimension 
t = - L ^ 9e ,~ 9z ) J 1//2 ^(0) — 2 and a dimensionless parameter y = where k is a cor- 

i- [ ?+(l9e|»?++9 z ) J w «(0) ' 

relation function given in ([21$, a = Jo°° ^r 3 n(r). Both t and y depend on the densities, 
through £ + = \{^- + and A = J^ll^p£_ ? respectively. Both the excitation gap in the 
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strong coupling regime and the density- dependent phase transition can be observed in ex- 
periments. On the theoretical side, it is interesting to make further studies of the model in 
the framework of bosonization [2, 
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